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Abstract 



For the first time the persistent current in a 2D free-electron system has 
been calculated analytically. The tight binding model is considered on a 
square lattice with filling factor 1/2. The array has a shape of rectangle with 
boundary conditions in both directions twisted by 2ircp x and 27n/> 3/ . The com- 
ponents of the twist are associated with two components of the magnetic flux 
in torus geometry. An analytical expression is obtained for the energy and for 
the components of the persistent current (PC) at a given flux and tempera- 
ture. It is shown that at zero temperature the PC density is proportional to 
the vector potential with the coefficient which does not depend on the size of 
the system. This happens because the Fermi surface for a square lattice at 
filling factor 1/2 is flat. Both the energy and the PC are periodic functions 
of the two flux components with the periods 4>o/q an d <fio/s where <po = hc/e, 
and q and s are integers which depend on the aspect ratio of the rectangle. 
The magnitude of PC is the same as in superconductors. Therefore, a 3D 
system constructed from a macroscopic number of isolated coaxial cylinders 
at zero temperature reminds the London's superconductor. It exhibits the 
quantization of trapped flux as well as the Meissner effect. However, all the 
phenomena are of a mesoscopic nature. The critical field H c decays with an 
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effective size of the system, H c ~ l/R e f- The magnitude of PC decays with 
T as ex.-p(—irTR e f/2at), where t is the hopping amplitude and a is the lattice 
constant. 

73.23.-b, 04.20.Jb, 71.45.-d, 75.20.-g 
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I. INTRODUCTION 



Persistent current (PC) in mesoscopic structures™!! has been extensively studied during 



between electrons™. 

The PC is a reaction of a system to an applied flux $, or, equivalently, it can be de- 
scribed as a change of the energy of the system due to twisted boundary conditions. In a 
two-dimensional system which forms a cylinder the twisted conditions mean that the wave 
function of a system acquires a factor exp(227r$/0o) with a circulation of one electron around 
the axis of the cylinder. Here 0o — hc/e is the flux quanta. 

The flux $ is related to the tangential component A of the vector potential $ = 2rrRA, 
where R is the radius of the cylinder. For a system with Galilean invariance the following 
simple statement is correct. The energy of a state with a given value of tangential component 
P of the total momentum depends on A as 



where N is the number of electrons and M = Nm is their total mass, m being the mass of 
one electron. The 2D current density j for a state with fixed P is 



where S is the area of the cylinder surface and n = N/S. At P = Eq. (|2]) reminds the 
London equation for a superconducting current. In this case n should be the superfluid 
density. 

A general derivation of Eq. (0), given above, is misleading because PC should be defined 
as a current in the ground state rather than in a state with fixed P. In two- or three- 
dimensional systems of free electrons the derivative of the energy with respect to A cannot 
be taken in such a simple way because the intervals of $, where branches of spectrum with 
different P change each other in the ground state, tend to zero with increasing system size. 





(1) 




(2) 
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For electrons in a periodic potential the situation is typically similar. The derivative of 
energy with respect to the flux is large for a given branch. However, different branches replace 
each other in the ground state at such small intervals of $ that the derivative taken at a given 
total quasimomentum P does not reflect properties of the ground state. Scalapino et al.0 
considered a tight binding model on a two-dimensional square lattice. Their computations 
show that at filling factor u = 1/4 the first level crossing occurs at $ ~ 1/L, where L is the 
size of the system. Their general conclusion is that the superfluid density, as found from the 
relation between j and A, is zero for free electrons in the tight binding model. We show in 
this paper that this is not always the case. 

Namely, we consider a 2D system of free electrons on a square lattice in a tight binding 
approximation at filling factor v = 1/2. The shape of the system is assumed to be a rectangle 
with arbitrary aspect ratio. We demonstrate below that at T = the two-dimensional PC 
density does not depend on the size of the system and has a form: 

j = ~—{A-A ). (3) 

Here m = h 2 /2ta 2 is the electron mass, t being the nearest-neighbor hopping energy. The 
two-dimensional density is determined as n = l/2a 2 , where a is the lattice constant. For 
simplicity, we consider a system of spinless fermions. The generalization to the case of 
non-interacting fermions with spin is straightforward. 

The constant A shows that the minimum of energy occurs not at zero flux. In contrast to 
Eq. (0), Eq. (|) describes PC in the ground state of the system which is a periodic function of 
$ with period 4>o/q. Eq. (|3]) is valid within the interval < § < (f) /q, or < A < <fi /27rRq, 
and is to be repeated periodically for other values of flux. Here q is an integer which depends 
on the aspect ratio of the cylinder and on the type of the boundary conditions imposed in 
the direction of the cylinder axis. The first term is shown to be independent on the aspect 
ratio. 

We found PC to have an order of magnitude of the London current. Note that this result 
gives substantially larger PC than it is prescribed in the ballistic regime by the so-called 



M-channel approximation (see Ref. |10| and references therein). Namely, our exact solution 
gives PC larger by a factor of \[L for the L x L square. This discrepancy is due to the fact 
that in the case of flat Fermi surface all transverse channels are coherent. 

Considering the 3D system constructed of a large number of coaxial closely packed 2D 
cylinders we show that it mimics the Meissner effect and the quantization of flux trapped 
in the opening. 

These properties appear since the Fermi surface at v — 1/2 is flat and no branch crossings 
occur in large intervals of $. Say, for a square array no branch crossing occurs in the whole 
interval < $ < <po, which means that q = 1. 

In fact, we are discussing a mesoscopic effect. The expression Eq. (Q) is valid only at 
mesoscopically small temperature, and the ideal diamagnetism occurs for mesoscopically 
small values of magnetic field: 

T < T c ~ -£-t, (4) 
a t 

where h is the spacing between neighboring coaxial cylinders. The effective size R e f is given 
by 



R ef = S D = q2nR = ^sq2irRD, (6) 

where D is the length of the cylinders and s and q are integers determined by the aspect ratio 
2nR/D (see below). Thus, this system can be classified as a "mesoscopic superconductor." 

Note that the average distance between energy levels in a 2D system is proportional to 
1/R 2 . The 1/R behavior in the above equations is also a result of the flat Fermi surface at 
v = 1/2. As it is seen from the calculations, all relevant interlevel distances are of the order 
of (a/R)t, rather than (a/R) 2 t. 

Since both T c and H c vanish at large R, there are no real critical phenomena in this 
model system. 
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II. CALCULATION OF PC AT ZERO TEMPERATURE 



Consider a rectangle of L x x L y lattice sites with periodic boundary conditions twisted 
in both directions by 27r$> x /(f)o and 27r$ y /0o- This corresponds to a toroidal geometry where 
$> x is the flux through the crossection of the torus and <& y is the flux through the opening. 

The single-electron energies have the form 

"2tt , 



e(n x , n y , (p x , <fi y ) = -2 J < cos 



27T 
Lin 



J X 



+ COS 



-r\ n v 



(7) 



where we introduce dimensionless <f> Xjy = Q X) y/4>o to simplify the notation. The values of 
integer quantum numbers n x and n y are restricted to the rectangle \n XtV \ < L x>y /2 (first 
Brillouin zone). To find the energy of the ground state one has to sum e(n x ,n y ) over the 
values {(n x ,n y )} inside the Fermi surface. 

In principle, the calculation of PC can be performed either at constant number of particles 
N or at constant value of chemical potential [l. Generally speaking, these two definitions 
are not equivalent. It is important to note that such a problem does not exist at v — 1/2 
at even numbers of L x ,L y at least. As one can see from the Eq. (0), every single-electron 
energy changes sign under the transformation n x , n y — > n x + L x /2, n y + L y /2. It follows, that 
at v — 1/2 due to the electron-hole symmetry the chemical potential \i is zero at any value 
of flux and at any temperature. Thus, if the flux changes at \x = 0, the number of particles 
in the ground state of the system does not change and if the flux changes at a given number 
of particles such that v — 1/2, the chemical potential does not change. 

Let us define the Fermi "surface" (FS) in two-dimensional n x ,n y space by the equation 

e(n x ,n y ,(f) x ,<f)y) = 0, (8) 

considering n x ,n y as continuous variables. It is easy to see that the FS forms a rhomb at 
any value of flux. Change in the flux produces a shift of the FS as a whole without changing 
its shape. 

First, let us consider for simplicity a square sample, L x = L y . The FS forms a square 
shown in Fig. 1(a). At <p x ^ y = some of the allowed single-electron states lie exactly at the 



sides of this square. All the states inside the square and 1/2 of the states at the sides of 
the square are occupied. All the states at the sides have the same energy so the occupation 
numbers of these states are not defined, while the many-electron ground state is degenerate. 

The degeneracy is lifted at infinitazimally small values of 0. Suppose that <p y = and 
<p x > 0. Then FS is shifted to the right (See Fig. 1(a)). All occupation numbers become 
defined. Namely, the states at the right side of initial square get occupied and those at the 
left side become empty. Note that the occupation numbers of as many as 2L states change 
when <p x crosses zero. 

It is easy to see that the occupation numbers are constant throughout the interval < 
4> x < 1. The total energy decreases with (ft x and then increases again. At <f> x = 1 all electrons 
jump one step to the right and the Fermi surface restores its original position with respect 
to the lattice of integer numbers (n x ,n y ). The total energy thus returns to the same value 
as at cj) x = 0. 

It follows that the total quasimomentum of electron system in the ground state does 
not change through all this interval and no branch crossing occurs. Then the sum over the 
occupied states can be easily evaluated: 

L/2-1 L/2-n y L/2 

E((f) x ,<f)y) = J2 E [e(n x ,n y ) + e(n x ,-n y )]+ e ( n x,0) (9) 

n y =1 n x =-L/2+n y +l n x =-L/2+l 
2m/L 



8t Re- 



e 



^ e 2m/L _ iy 

This expression is exact in the region < (f) x ± (j) y < 1. In the limit of large L, the 0- 
dependent part of the energy, SE((p x , <p y ) = E((f> x , (f) y ) — E(0, 0), can be written in the form: 

8E{<t> x , <j) y ) = 8* Uj - 0.(1 - (j> x )\ . (10) 



Repeating Eq. (|T^) periodically one gets the expression valid in the whole plane 



6E(<f> x , <j)y) = At 



1\ 2 , 1\ 2 1 



^-2j + 1^-27 2 



(11) 



where 0± = (f> x ± (fi y , and {...} denotes the fractional part of the number, defined as a 
difference between the number and the largest integer less than it. 



Fig. 2 shows the energy SE(<p x , 4> y ) as given by Eq. (|TT|). The positions of energy minima 
form a square lattice shifted from the origin: 

i-j 



(12) 



2 2 
with arbitrary integer % and j. 

The point <p x = <p y = corresponds to a maximum of energy, in the same way as 
in ID case with even number of electrons. At this point the derivatives d5E / 'd(p XyV are 
discontinuous. Such behavior appears as a result of lifting of the degeneracy of the states 
at the Fermi surface. 

The PC at T = can be found as the derivative of the total energy with respect to flux: 



dE 



<9$ 

00 



x,y 



_c_ f dE dE \ 

0o \^0+ 5 0-/ 



(13) 



The magnitude of 6E((j) x ,<f> y ) and I x>y ((f> x ,4> y ) as given by Eqs. fly} ), ( |l3l) is independent 
on the size L of the square. Such a large magnitude results from the fact that in the region 
with no branch crossings (or, with no electron changing its state) all electrons together 
contribute to the current. 

It may seem that the aspect ratio L x /L y = 1 is crucial for the effect. In the next section 
we calculate PC at finite temperature for arbitrary aspect ratio L x = sK, L y = qK with 
mutually-prime integers s and q. We assume macroscopic limit K — > oo. It is useful to 
generalize <f>± for a rectangular sample as 



± S( 



(14) 



In the limit T = we find 



8 ct 
sq 0o /< 



8 ct 
sq (fio/s 
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{0 + } 
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{0-} 



(15) 



5E((j) x ,(j)y) ~ 

sq 



(16) 



The flux- dependent part of the energy can be restored from Eq. (jig): 

^ + >-i) 2 + (<*->4) 2 -[ 

This result is a generalization of Eq. (JTTJ) to an arbitrary aspect ratio s/q of the rectangular 
sample. 

As follows from Eqs. (|1|), flTB]), the energy and current as functions of flux do depend 
on the aspect ratio. However, they do not depend on the system size, if the aspect ratio is 
kept constant. 

The result Eqs. (|l"5l), fllTf), can be understood from Fig. 1(b), which is drawn for the 
case L x /L y = 2/3. Contrary to the Fig. 1(a), there are now points (n x ,n y ) closer to the 
Fermi surface than one lattice spacing. However, there is still a regularity in their positions. 
Namely, as the Fermi surface shifts with flux, the points enter the Fermi sea in groups. 
Consider, for example, the same case as above: <p y = and (p x > 0. As seen in Fig. 1(b), the 
branch crossings occur only at (f) x = 2ir/3, 47r/3, and 27r. In terms of <fi± this corresponds 
to integer <f>± = 3(f) x = 1, 2, and 3. These values of flux are determined by s and q and 
do not change with the size of the system. The number of points in each group, in tern, is 
proportional to the size of the system, so the corresponding contribution to the current is 
large. 

At $ y = Eq. fllD gives 

Ix = ~~qWv{{j£\~^' (17) 
Up to now we are discussing the torus geometry. To come to a cylinder geometry one 

has to formulate the boundary conditions in the direction of the cylinder axis, chosen as y. 

In what follows we assume periodic boundary conditions in this direction with § y = 0. This 

leads to Eq. ([Tj]) for a total current through the cylinder. As another option me may impose 

the condition that the wave function is zero at the edges of the cylinder. It can be shown 

that in this case the second term in Eq. ( |TTD changes while the first term remains intact. 

Note that both energy and current are periodic functions of flux with period <po/q rather 

than 0o- 
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Taking into account that the current density j x = l x j (aL y ) and that the vector potential 
Ac — &x/( a L x ) one obtains Eq. with the first term independent of s and q. 



III. PC AT FINITE TEMPERATURE 



We start with the equation 

j_ c^^ de(n x ,n y ) 1 

<h£o£o <90, l + exp(e(n„ % )/T)' 

It is convenient to rewrite the single electron energy in the form 

e(n x , n y ) = -At cos {j^( n + ~ <t>+) \ cos {j^{ n - ~ <t>-)^ > ( 19 ) 

where n± = gn, ± sn y and 0± are given by Eq. (|H|). Using d/d(p x = q(d/d(p + + <9/<90_) we 
find that the current has two terms, 

I x = - s (I + + I-), (20) 



where 



The idea of our calculation is to transform Eq. fl2"T| ) in such a way that the internal sum 
gives PC of ID problem with effective temperature and effective flux. For this purpose we 
use the identity: 

sK-lqK-l sK-1 q-1 K-l 

E E /("». ny) = E E E /("» + rf + **, d + qk). (22) 

n x =0 n y =0 m=0 d=0 k=0 

This identity is valid for any function f(n x , n y ) periodic in n x and n y with periods sK and 
qK respectively. Then I + can be written in the form 

Airct s ^ ^ K- 1 sin f-^(n + - 0+)) cos (-^(n. 



E E E r 4 / V£9£J_ — ^ ry. (23) 

"vd m=o d=o k=o 1 + exp [— cos (^(n+ - 0+)J cos (^(n_ - 0_)J 
Similar expression can be written for 7_. Note that n_ = gm + (q — s)d does not depend on 
k, while n + = (qm + qd + sd) + 2sqk does depend on k. Therefore, the current I + can be 
written as 
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sK-1 q-1 

i+= E ( 24 ) 

m=0 d=0 

where I + (m, d) denote the internal sum over k, 

if*, T ^ 1 + exp [-f cos (f (k - + ))] 

The sum in Eq. ( p5|) describes the PC in ID system with effective temperature and effective 
flux given by 

T 



T(ra, d) = -, 



<p+{m,d) 



(26) 

qm — qd — sd 



2sq 



Using the Poisson summation formula (see Ref. |15[), one obtains 



A 8ircT ~ cos(IitK/2) . . ~ . 



Performing the summation over m and d in Eq. (|24[ ) we note that T is a smooth function 
of m/K and d/K. However, sin(27r/0 + ) has an oscillatory behavior for some I, so that the 
contribution of the corresponding harmonics vanishes in the limit K — > oo. The oscillatory 
behavior is absent if I is an integer multiple of 2sq. For these I, the sum over m can be 
transformed into integral via p = (7r/sK)m, while the sum over d simply gives a factor q. 
Thus, one obtains 

oo 

I ± = X>;sin(27d0 ± ), (28) 
i=i 

where 

A l = sqKT^ f . . (29) 

(po Jo smh{L7rsqK 1 / 2tsmp) 

For the PC in x-direction one has from Eq. ( p0|) 

J oo 

I x = A i [sin(27r/0+) + sin(2vr/0_)] . (30) 
s i=i 
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Similar calculation gives 

oo 

I y =-^TAi [sin(27r/0 + ) - sm(2vr/0_)] . (31) 
9 1=1 

Eqs. (|2TS|), ( POD give the Fourier series expansion of PC at any temperature. Expansion 
of Ai at small KT/t yelds 

A, « £1 (32) 
In this case Eqs. (p0|), (|31|), and ( p2|) give the Fourier series expansion of the zero-temperature 



result Eq. (|T5|). 



In the opposite limit, KT/t 3> 1, the amplitudes of the harmonics decay as 



lirsqKT\ 



sqKTtexpl ^ J , (33) 



00 V 7 / 

so that PC is dominated by its lowest harmonic. When <& y = one has 

sq 0o/ 



6c ( R ef Tt \ 1/2 ( nR ef T \ . f ® x \ 



IV. LOW TEMPERATURE MAGNETIC PROPERTIES 

In this section we study magnetic properties of a quasi-3D system constructed of a 
macroscopic number of closely packed coaxial cylinders assuming that the temperature is 
very low. Then the connection between flux and current for each cylinder is given by Eq. (|17|) . 
For the sake of simplicity we assume that the cylinders are long, such that the circumference 
of the internal cylinder 2ttR = aL x is much larger than D = aL y . The distance between the 
internal and external cylinders is supposed to be much less than R. We assume further that 
all the cylinders have the same ratio L x /L y = s/q. One can imagine a small change either 
in L x and L y of adjacent cylinders or in their lattice constant. 

The second term in Eq. (fi7| ) appears since zero flux does not correspond to the minimum 
of energy. It may lead to an appearance of a spontaneous flux in this system. This idea has 
been put forward by Wohlleben et al, and Szopa and Zipper, Ref. [l^, and then studied in 
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details in Ref. |13| These authors considered a cylinder constructed from isolated ID rings. 
Loss and Martin0 argued that in a single ID ring no symmetry breaking can occur, but 
their arguments are restricted to ID case. 

In this paper we concentrate on the first term in Eq. fllTD - It is an analog of the Lon- 
don current in superconductors and it creates a strong diamagnetism in a quasi-3D system 
described above. Suppose that an external magnetic field H ext is applied to the system and 
that there is a solenoid creating flux $> ext inside the internal cylinder. 

Let $fc be the total flux inside cylinder k, where k — 1 for the internal cylinder and 
k = N for the external one. The flux obeys the equation 

/ 4tt n \ 

$ fc - = 2nRb H ext + — £ . (35) 

Here b is the distance between adjacent cylinders which we assume to be of the order of the 

lattice constant a. Since the thickness d = Nb is supposed to be much less than R we have 

neglected that the radii of cylinders are slightly different. The right hand side of Eq. (|35D 

describes the flux through the area between the k-th and (k — l)-th cylinders created by the 

external field and outer cylinders. 

The following condition should be added to this finite difference equation: 

/ 4vr n \ 

$! - <$> ext = ttR 2 \H ext + m) . (36) 

If $A: is a smooth function of k one can transform Eq. (|3^) into differential equation 

rf 2 $ (j) /q (\ $ 1 1\ 



dr 2 " X 2 \\ (p /q J 2/ 
Here A is the analog of the London penetration depth 



(37) 



2 4vrl6t 16e 2 n 3 



b (J)q Time 2 
where = l/2ba 2 is the 3D electron density. 

Eq. (|36| ) transforms into the boundary condition at r = R: 
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(39) 

--R 



The second boundary condition reads 

= 2nRH ext . (40) 



dr 



r=R+d 

One can use differential equation if A b. 

Eq. ( JITD can also be obtained by minimizing total energy with respect to flux. The total 
energy consists of two parts. First is the energy of magnetic field in the space between 
cylinders. The magnetic field can be expressed through d§/dr using Eq. ( |35l) as 

^ = 2nRH(r). (41) 
dr 

The second part is the internal energy of 2D electron gas. This energy per cylinder is 
given by Eq. (|1^) at <p y = 0. Thus, one gets for the total energy 



E, 



1 D /(§)' w> 



/ ol III 



8?r 2-kR 

Minimizing this expression with respect to $(r) and taking into account that d5E/d& = 
(— l/c)/($), where /($) is given by Eq. (|TTD, one obtains Eq. (|T7|). 

The Eq. (|37|) is nonlinear since it contains the fractional part {$/(0o/<?)} which makes 
the right-hand side periodic. However, it becomes linear if the total drop of the flux inside 
the system is smaller than <fio/q. If H ext = the solution of the linearized equation with 
boundary conditions (|39|) , (f40D in the case R 3> d ^> A is 

$(r) = $ n + {§ ext - $„)^ exp f- 1 -^) . (43) 



R r V A 
Here 

*« = 7("-5)- (44) 

One can see that the flux inside the cylinder with d 3> A may take only quantized values 
$ n with arbitrary integer n. Note that there is no zero flux among the allowed values of 
the frozen flux $„. This is because zero flux does not correspond to a minimum of the total 
energy at zero temperature. The solution Eq. (f43D is obtained in the linear approximation 
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and it is valid if ($ ea;i — $ n )2A/i? < 4>o/q. The physics of this result is that the inner cylinders 
carry a current which creates a favorable flux for the rest of the system. 
If the system is in an external magnetic field H ext , the solution is 



In this case the cylinders near external surface carry current which screens magnetic field 
inside the system and adjusts the total flux to The solution is valid if 2nR\H ext < 4>o/q- 
This condition is equivalent to Eq. (|5]). It has a simple interpretation. The loss in the total 
energy due to the ideal Meissner effect is of the order of H^ xt RDb per cylinder. The gain in 
the energy of a cylinder due to the adjusted flux is of the order of t/sq (see Eq. (|16D). At 
large field the loss becomes larger than the gain and the field penetrates into the system. 
This is the origin of a "mesoscopic" critical field. Note that the relation H^RDb ~ t/sq is 
also equivalent to Eq. (|5|). 

It follows from the results of the previous section that zero-temperature approximation 
is good if sqKT/t = Ty/sq27cRD/at <C 1. This is the same condition as Eq. (|j). At larger 
temperatures the penetration depth A increases as exp(irTR e f /4at) and eventually reaches 
the thickness d of the cylinder, gradually destroying strong diamagnetism. 



Finally we have presented a model which mimics in a mesoscopic scale some properties 
of superconductors, such as Meissner effect and quantization of flux, though the physics of 
the model does not involve any electron pairing. The flux quanta in the model is <fio/q where 
q is determined by the aspect ratio of the system. 

Since the range of temperature and magnetic field for these phenomena shrinks to zero 
in a macroscopic system, one should not expect any phase transitions. However, for a 




(45) 



or in terms of magnetic field defined by Eq. ([41]) 




(46) 



V. CONCLUSIONS 
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mesoscopic system this range is not necessarily small. Let us assume a hypotetic 3D layered 
system with very weak interaction between layers and flat two-dimensional Fermi surface. 
Then it follows from Eqs. (^), (|5]) that the temperature range is up to 12K and the range 
of H ext is up to 240 gauss for a system with R e f = 3 x 10 _5 cm, a = b = 3 x 10~ 8 cm, and 
t = leV. In a system with disorder the obvious condition for these phenomena is that the 
elastic mean free path is smaller than the size R e f. 

Our model ignores electron-electron interaction. We hope that it is not important at 
large t. Our modeling of small interacting systems up to 18 electrons show the same value 
of the PC at t immediately above the Wigner crystal quantum melting point0. 
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FIGURES 

FIG. 1. The Fermi surface at v = 1/2. The points represent allowed integer values of n x and 
n y inside the first Brillouin zone. The dashed lines show the Fermi surface at zero flux. The solid 
lines are the Fermi surface shifted by flux. The aspect ratio L x /L y = 1 (a) and L x /L y = 2/3 (b). 

FIG. 2. Lines of constant SE(cj) x , <f> y ) as given by Eq. ([[(]) for square sample s = q = 1 (a), and 
for rectangle with s = 2 and q = 3 (b). Note the difference in periodicities. 
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